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ABSTRACT 


This  thesis  is  concerned  with  the  steady  state  response  to  a  step 
input  of  a  non-linear  servomechanism  having  viscous  friction  and  a  load 
position  unity  feedback  loop.  A  gear  train  with  backlash  between  re¬ 
silient  gear  teeth  was  located  between  the  motor  and  load. 

The  physical  equations  describing  the  system  are  presented  and 
adapted  to  a  form  for  the  Control  Data  Corporation  l6o4  digital  com¬ 
puter.  The  computer  program  used  in  this  analysis  is  presented  along 
with  its  flow  diagram  and  description  of  its  operation. 

The  results  of  this  thesis  are  presented  in  several  forms  using 
the  parameters  of  system  damping  coefficient,  distribution  of  friction 
and  inertia,  backlash  angle  and  coefficient  of  restitution.  The  re¬ 
sults  are  applied  to  sample  problems  in  illustration  of  their  use  for 
design  and  analysis. 

The  authors  wish  to  express  their  appreciation  to  Dr.  George  J. 
Thaler  of  the  Department  of  Electrical  Engineering,  and  to  Dr.  William 
Wainwright  and  Mr.  Edward  N.  Ward  of  the  Department  of  Mathematics  for 
their  assistance  and  encouragement  in  completing  this  work. 
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*(  )  Terms  indicate  computer  mnemonics  referred  to  in  Sec.  3,  Computer 
Program  Development. 
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1.  Introduction  and  Background 

The  problem  under  investigation  was  the  response  to  a  unit  step 
input  of  a  nonlinear  servomechanism  having  viscous  friction  and  load- 
position  unity  feedback  loop.  A  gear  train  with  backlash  between  the 
resilient  gear  teeth  was  located  between  the  motor  and  load.  A  de¬ 
termination  of  the  possible  existence  of  a  limit  cycle  with  variation 
of  several  parameters  was  made.  The  variable  parameters  considered 
were  motor ,  load,  and  system  time  constants,  backlash  and  the  co¬ 
efficient  of  restitution  of  the  gear  teeth.  If  a  limit  cycle  existed, 
a  further  investigation  wad  made  into  the  change  of  the  size  of  the 
limit  cycle  with  variation  of  the  above  parameters. 

If  a  gear  train  is  required  between  the  motor  and  load  in  a 
servomechanism,  the  gear  train  may  be  treated  as  ideal  with  perfect 
meshing  of  the  gears.  In  this  case  the  components  of  the  system  are 
joined  at  all  times  and  the  system  may  be  described  by  a  single  linear 
differential  equation.  The  classical  linear  solution  in  either  the 
frequency  or  time  domain  is  the  result. 

Present  production  methods  cannot  meet  the  requirements  of  an  ideal 
gear  train,  having  perfect  engagement  of  the  teeth,  nor  is  the  ideal 
gear  train  desirable  from  the  standpoint  of  wear  on  the  gear  faces. 

The  result  is  that  the  practical  gear  train  has  separation  of  the  gear 
teeth  when  the  velocities  of  the  motor  and  load  are  different.  The  re¬ 
sponse  of  a  servomechanism  with  backlash  is  in  effect  discontinuous  or 
nonlinear.  The  net  system  acts  as  three  individual  but  dependent  systems 
with  the  boundary  conditions  of  each  being  specified  by  certain  physical 
laws . 

Two  accepted  methods  of  analysis  are  available  for  solution  of  the 
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problem.  The  first  is  the  solution  by  describing  function  methods 
while  the  second  is  by  analytical  methods.  Describing  function 
methods  are  primarily  frequency  response  techniques  and  depend  upon 
the  assumption  that  the  nonlinear  element  may  be  considered  as  linear 
over  the  range  of  consideration.  A  further  assumption  is  made  that 
the  input  to  the  nonlinear  element  is  a  pure  sinusoid  with  no  harmonic 
frequencies.  The  describing  function  method  is  primarily  concerned 
with  steady  state  results. 

The  analytical  solution,  which  is  the  second  method  of  considering 
the  problem  is  primarily  a  ^transient  response  technique  concerned  with 
solution  in  the  time  domain.  This  approach  also  requires  certain 
assumptions  concerning  the  behavior  of  the  system.  The  assumptions 
used  in  several  investigations  are  pointed  out  in  the  following  para¬ 
graphs,  while  the  assumptions  of  this  investigation  are  outlined  in 
Section  3>  Development  of  Equations.  The  validity  of  the  analysis  by 
any  method  depends  on  the  accuracy  with  which  the  real  system  has  been 
defined. 

Chestnut  and  Mayer,  Ref.  a,  have  considered  the  backlash  problem 
using  describing  functions  and  supported  the  findings  by  analog  com¬ 
puter  studies.  The  main  difference  between  this  cited  work  and  previous 
describing  function  analyses  was  the  simulation  of  backlash  by  a  dead 
zone  between  the  motor  and  load  and  springiness  in  the  interconnecting 
shaft . 

Various  analytical  approaches  have  been  made  to  the  problem  of 
backlash,  Lutkenhouse,  Ref.  b,  used  graphical  methods  in  the  phase 
plane  to  solve  several  cases  of  plastic  impact  between  the  gear  teeth 
of  a  second  order  servomechanism.  Pastel  and  Thaler,  Ref.  c,  developed 
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analytical  equations  for  the  existence  of  limit  cycles  using  phase 
plane  equations  for  the  case  of  plastic  impact  and  no  load  inertia. 

Knoll  and  Narud,  Ref.  d,  investigated  limit  cycles  in  the  phase  plane 
using  an  analog  computer.  The  investigation  of  Knoll  and  Narud 
covered  a  wide  range  of  parameters  for  the  case  of  plastic  contact 
between  the  gear  teeth. 

To  analytically  describe  the  total  system  response  of  a  second 
order  servo  with  backlash,  three  differential  equations  and  one  or 
more  algebraic  equations  are  required.  The  system  as  a  whole  may  be 
treated  as  being  piece  wise  linear.  One  differential  equation  is  used 
to  describe  the  entire  system  when  the  gear  teeth  are  in  contact  and 
the  motor  is  driving  or  braking.  When  the  gear  teeth  are  not  in  con¬ 
tact,  two  more  differential  equations  are  required,  one  for  the  motor 
alone  and  one  for  the  load  alone.  The  boundary  conditions  of  the  dif¬ 
ferential  equations  are  determined  from  the  solution  of  one  or  more 
algebraic  equations,  expressing  the  laws  of  conservation  of  momentum 
and  energy.  For  plastic  impact  with  no  bounce  of  the  gear  teeth,  the 
law  of  conservation  of  momentum  must  be  satisfied.  For  perfect  elastic 
contact  between  the  gear  teeth  the  law  of  conservation  of  momentum  and 
the  law  of  conservation  of  energy  must  be  satisfied  simultaniously. 
Intermediate  cases  between  perfect  elastic  and  perfect  plastic  contact 
can  fulfill  only  the  law  of  conservation  of  momentum.  However,  total 
energy  accounting  may  be  made  in  the  intermediate  and  perfect  plastic 
impact  cases. 

New,  Ref.  e,  adapted  the  differential  equations  of  system  motion 
and  the  law  of  conservation  of  momentum  (plastic  contact)  to  a  digital 
computer  analysis  of  a  second  order  servomechanism  with  backlash.  The 
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digital  computer  was  chosen  as  the  method  of  solution  for  this  investi¬ 
gation  primarily  "because  of  the  flexibility  of  data  presentation  and  the 
ease  with  which  cyclic  parameter  variation  could  be  obtained. 

In  later  sections  of  this  thesis,  the  differential  equations  for 
the  system  and  the  algebraic  equations  for  the  impact  boundary  conditions 
are  developed.  The  computer  program  and  its  flow  diagram  and  modes  of 
operation  are  pointed  out.  The  results  of  the  solutions  of  problems 
are  presented  in  several  forms  and  observations  are  made  in  Section 
Results  and  Discussion.  Sample  applications  of  the  results  to  mechani¬ 
cal  systems  are  presented  fn  Section  6,  Application  of  Results. 

At  this  writing,  associated  work  in  the  area  of  transient  response, 
primarily  peak  o-.  ershoot  and  settling  time,  of  a  second  order  servo¬ 
mechanism  with  backlash  and  plastic  and  elastic  impact  is  being  pre¬ 
pared  as  a  thesis  by  C.  E.  Andrews  and  R.  A.  Kelley  at  the  U.  S.  Naval 
Postgraduate  School. 
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2.  Development  of  Equations 

The  equations  developed  were  those  of  a  phase  plane  analysis. 

Prior  to  defining  the  net  system  equations,  the  assumptions  under  which 
the  analysis  was  made  will  be  stated.  When  necessary  these  assumptions 
will  be  amplified  and  referred  to  later  in  this  work. 

It  was  assumed  that: 

1.  The  gear  teeth  were  initially  in  contact  and  the  initial 
conditions  of  the  system  were  all  equal  to  zero.  This  was 
later  proved  to  be  an  unnecessary  limitation  for  the  study 
of  steady  state  response. 

2.  Plastic  deformation  of  the  gear  teeth  during  steady  contact 
and  impact  and  any  torsional  deformations  of  driving  shafts 
are  negligible. 

5.  The  inertias  of  the  gears  and  drive  shafts  are  considered 
as  part  of  the  load  or  motor  inertia  depending  on  their 
attachment  in  the  system. 

4.  The  law  of  conservation  of  energy  was  completely  satisfied 
in  only  the  perfect  elastic  case  by  maintaining  the  total 
mechanical  rotational  energy  of  the  system  constant  at  the 
instant  prior  to  and  after  impact.  The  law  of  conservation 

of  momentum  is  satisfied  in  plastic,  elastic,  and  intermediate 
cases.  When  the  law  of  conservation  of  momentum  is  the  only 
equation  required  to  be  satisfied,  the  energy  lost  from  the 
system  is  dissipated  in  the  heat  of  infinitesimal  deformations 
of  the  gear  teeth. 

5.  The  gear  teeth  are  in  contact  only  instantaneously  during 
impact  for  the  elastic  and  intermediate  cases  (excluding 
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plastic)  and  that  the  impulse  torques  of  drive,  friction  and 
hearing  supports,  etc.,  are  zero  during  impact,  e.g. 

j  Tm  —  O  CXVAcl  Qf  —  0- 

6.  The  coefficient  of  restitution  of  the  two  opposing  gear 
teeth  is  the  same  or  is  described  by  an  equivalent  co¬ 
efficient  if  the  two  gear  teeth  are  of  unequal  coefficients. 

7.  Backlash  is  assumed  to  be  equal  at  all  points  on  the  gear 
circumference.  Backlash  is  measured  at  the  output  shaft. 

A  block  diagram  of  the' system  considered  is  presented  below. 
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^  =  -  p 
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and  factoring 


The  momentum  equation  which  is  independent  of  e  is  ;  *  ^ 

lift,  f  --  (6'c- 

p  ^  u  V  -  »  /  v  /  '  /  *  \ 

or  rearranging  ^  J  “  (  £>C~  j 

The  factored  energy  equation  is  now  divided  by  the  momentum 

equation  to  yield  ^  ^  ~  /  "&  £  ~h  &£,  ) 

which  checks  with  the  definition  of  the  coefficient  of  restitution 

for  e=l. 


The  momentum  equation  implied  from  the  impulse  approach  that 
dH:-*  0  hence  Qc=  &c  0\Y\h  .  The  definition  of  the 

coefficient  of  restitution  is  independent  of  position.  It  may  he 
noted  that  the  equations  for  energy  and  momentum  depend  only  on  the 
inertia  distribution,  gear  ratio  and  angular  velocities. 

In  manipulation  of  these  equations  it  is  seen  that  load  inertia 
may  be  transferred  to  the£|  ,  shaft  by  multiplication  ^  ,  the  same 

transfer  may  be  accomplished  with  load  friction.  Since  ^  ,  and  , 
are  related  to  and  it  may  be  reasoned  that  the  results 
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for  =  1  may  be  extrapolated  to  physical  systems  where  /=  1. 

Since  a  total  energy  accounting  may  be  made  for  the  system  with 
the  equations  used,  the  rotational  energy  lost  is  attributed  to  the 
heat  of  deformation  of  the  gear  teeth.  This  transfer  of  energy  could 
in  fact  be  determined  from  the  equations  used.  Thus  the  conservation 
of  rotational  energy  e=l  is  a  special  case  of  a  broad  interpretation 
of  the  law  of  conservation  of  energy. 
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The  equation  for  the  motor  alone  with  the  gears  not  in  contact  is. 

Tw  Qyv\  4*  f  yy)  ~  ~ 


The  method  of  analysis  used  by  New,  Ref.  e,  is  used  in  order  to 
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and  introducting  the  slope  equations  of  the  phase  plane 
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If  the  load  equation  is  first  put  in  the  form  for  the  phase  plane 
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The  same  result  may  be  obtained  by  setting  4),  =1.  The  results  of  such 
a  transformation  require  inverse  scaling  for  practical  application, 
examples  of  which  given  in  Section  6.  It  is  noted  that  the  slope  of 
the  load-free  equation  is  mathematically  the  same  whether  in  the  phase 
plane  or  in  the  transformed  phase  plane.  It  is  pointed  out  at  this 
time  that  the  equations  later  developed  to  satisfay  the  laws  of  con¬ 
servation  of  momentum  and  energy  are  independent  of  the  system  natural 
frequency. 

To  establish  equations  for  the  law  of  the  conservation  of  momentum 
and  a  relationship  satisfying  the  law  of  the  conservation  of  energy, 
the  following  schematic  is  used. 


It  is  seen  that  ^  ^ \  )  & C.  ~~  ^ I 

_/  , 

and  since  P  - 1  an  equivalent  schematic  results 


By  representing  the  inertia  of  the  load  and  motor  as  inertia  of  the 
gears  the  following  figure  is  obtained  at  the  instant  of  impact. 
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The  torque  equation  may  be  expressed  \  =  T9  -  Tdj> 

(X  t 

Impulse  is  equal  to  the  time  rate  of  change  of  momentum.  The  ex- 
pressions  for  the  rates  of  change  of  momentum  of  the  gears  treated 
separately  may  he  written: 
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where  the  primed  velocities  are  those  following  impact. 

If  the  momentum  is  to  he  conserved  in  the  system  the  impulse 
functions  are  equal  and  opposite. 

3J2L  ( M  ( K  -  &<) 

I  * 

and  3m  J_  (  -  6^)  -  (  6-c  -  ) 

P 

from  substitution  of  the  equations 

r.'  =  r»  ,  ©c  =  '(° 

The  assumption  that  all  other  impulse  functions  are  zero  at  the  instant 
of  impact  is  restated  at  this  point. 

A  definition  of  the  coefficient  of  restitution  will  he  made  for 
purposes  of  this  work: 

•e.  -  —  ( ■=.  —  (  P  6^  -  ^e, ) 
(4^  ■*-  4^ 

It  will  he  shown  that  this  definition  of  e  will  satisfy  the  law 
of  conservation  of  energy.  Consider  the  case  of  e=0,  f 

6--( pi'- 

in  which  case  for  plastic  impact  (e=0)  the  gears  are  moving  at  the 
same  velocity  following  impact. 

The  expression  for  the  coefficient  of  restitution  is  examined  for 
the  case  of  perfect  elastic  contact,  e=l,  for  which  the  law  of  conser¬ 
vation  of  energy  is  satisfied. 

Solution  of  the  restitution  equation  for  e=l  yields 
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3.  Computer  Program  Development 
The  physical  equations  of  the  net  systems  were  programmed  for 
solution  using  the  Control  Data  Corporation  1604  high  speed  digital 
computer  utilizing  paper  tape  program  input  and  magnetic  tape  output. 
An  IBM  717  line  printer  was  used  to  extract  data  from  the  magnetic 
tape.  The  Control  Data  Corporation  machine  library  and  the  U.  S. 

Naval  Postgraduate  School  computer  subroutine  library  were  used  for 
assembly,  Runge-Kutta-Gill  numerical  integration,  and  decimal  output. 

Several  changes  were  made  in  the  forms  of  the  physical  equations 
of  the  net  system  in  order  to  eliminate  duplicate  computing  operations 
and  to  fit  the  equations  to  a  form  suited  to  the  variable  parameters. 

In  the  table  of  symbols,  computer  mnemonic  (m  )  terms  which  are  used 
in  the  assembly  subroutine  have  been  indicated  by  parenthises  and  will 
be  defined  when  encountered. 

The  equation  of  the  system  with  motor  and  load  combined 

=  UV  -  £u 

was  put  in  the  form 

‘  *  7,  • 

-  COy\  ■Ac.  -_A  6c.  wkerc 

A  =  and 

THETAL  =  ©  r  ONE  = 

THETALD  =  U  -  ZETA  =? 

*  » 

UDOT  =  0^  QMEGANSQ 

were  the  computer  m  terms  used.  Underlining  will  be  used  in  this 
section  to  denote  the  transition  to  the  m  terms  of  the  computer 
equations. 

The  differential  equations  were  solved  using  the  Bunge-Kutta-Gill, 
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(RUNGE),  numerical  integration  method-  The  RUTJGE  subroutine  re¬ 
quired  the  definition  of  synonymous  terms  for  the  four  iterative 
cycles  used  to  produce  one  extrapolated  set  of  variables.  The  incre¬ 
ment  of  the  independent  time  variable  chosen  was  0.01  sec.  for  all 
0.1  and  0.004  sec.  for  all  ^  ^  6  •  |  • 

The  load  free  equation 


m  #  r 

-1-  £l  B  c  “  o  was  put  in  the  form 

*  *  Lm  * 

Be  4-  G  °  for  the  computer. 

A  separate  set  of  RUNGE  synonyms  was  used  for  the  load  free  equation 
in  order  not  to  destroy  the  previous  computations  before  they  were 
determined  of  no  further  value  in  computing  and  also  in  the  making 
of  program  decisions. 


(~)^was  designated  V  and  Q  was  designated  VDOT  for  the  computer 


program . 

The  equation  for  the  motor  alone 

Om  +  ^ w\  +  &c-  -  JS. 

T  ^  0  'Ad  1 

was  put  in  the  form 

^  (&R-  &c.) 

d  rYl 

-  -  b.  6  +  C  ( 6 £  - 


and  for  RUNGE,  the  terms 
©fh  =  THETAM 


YA 


=  THETAIM  =  W 


and 


6 


YY\ 


=  WDOT 


were  used. 


The  equation  representing  the  law  of  conservation  of  momentum  and 
the  definition  of  coefficient  of  restitution  were  combined  to  the  forms 


l6 


7  XtX  —  \  P  &  KVa  (  *  ^  ^  +  -t) 


•  I 


6. 


Vw 


for  the  computer  program.  No  provision  for  additional  m  terms  was 
made  to  denote  the  primed  values.  The  additional  terms: 


e  -  Rest  ituT  p  =  RHo  and  0*=  RHosq 


were  used. 

Two  equations  were  used  to  define  the  boundaries  of  operation  of 
the  load  and  motor  when  they  were  acting  separately, 


—  p  yy\  (XV\  d  0  p  &YV\  ^  A 


Two  major  decisions  of  the  computing  cycle  were  the  determination  of 
the  point  where  the  motor  and  load  would  float  free,  FLQATEST,  and 
the  response  of  the  system  to  impact  of  the  load  and  motor  gears  when 
the  boundary  conditions  on  position  were  met,  COMBTEST. 

The  first  decision,  the  point  of  float  free,  is  made  on  the 
equality  or  inequality  of  slopes  in  the  phase  plane,  i.e.,  at  separation 
Nl=  Ng.  Equating  the  slopes: 


t  B  =■  6 


The  above  equation  is  solved  after  computing  each  point  in  the 


combined  phase  trajectory.  If 


()  the  system  re¬ 


mains  combined.  If 


having 


the  response  motor  alone  and  load  alone  immediately  thereafter. 


The  second  major  decision,  that  of  whether  after  impact  the 
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response  should  he  that  of  the  combined  system  or  the  response  of  the 
load  and  motor  acting  separately,  is  made  on  the  basis  of  the  resultant 
velocities  after  impact  and  upon  which  side  of  the  backlash  the  motor 
and  load  positions  are  found. 

The  computer  program,  its  flow  diagram,  a  description  of  the  pro¬ 
gram  operation,  and  an  explanation  of  the  readout  of  information  with 
examples  are  presented  in  the  Appendix. 
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4.  Methods  and  Scope  of  Examination 

Two  modes  of  operation  of  the  computer  program  were  used.  The 
first  mode  considered  was  the  phase  trajectory.  This  mode  was  used 
to  check  the  operation  of  the  computing  cycle.  Time,  velocity  and 
position  of  the  load  were  printed  when  the  system  was  comhined.  When 
the  system  was  separated,  the  values  of  motor  velocity  and  motor 
position  were  also  printed.  Although  the  time  increment  used  for 
computations  was  either  0.01  sec.  or  0.004  sec.,  printed  outputs 
for  the  first  mode  were  taken  only  every  0.1  sec.  Additional  printed 
outputs  were  taken  at  the  time  of  contact  of  the  gears  just  prior  to 
and  immediately  following  the  solution  of  the  momentum  and  restitution 
equations.  A  typical  printout  of  the  phase  trajectory  mode  of  operation 
is  shown  in  Appendix  C. 

The  second  mode  of  operation  was  used  for  printing  only  the  maxi¬ 
mum  computed  positive  load  overshoot  position  for  each  cycle,  the 
associated  load  velocity  and  the  exact  problem  time  of  the  computation. 

A  sufficient  number  of  print-outs  for  each  problem  was  obtained  to 
prove  limit  cycle  existence  and  average  size.  Since  this  investigation 
was  concerned  primarily  with  the  steady  state  response,  the  second  mode 
of  operation  was  the  one  utilized  to  obtain  the  majority  of  the  data. 
This  read-out  method  markedly  decreased  the  data  reduction  time  for  the 
problem  analysis,  as  opposed  to  the  analysis  of  steady  state  conditions 
provided  by  a  full  phase  trajectory  print-out.  Since  the  time  required 
for  computer  read-out  was  several  magnitudes  greater  than  the  computing 
time,  valuable  computer  time  was  saved  by  this  mode  of  operation.  An 
example  of  the  printed  output  obtained  when  only  maximum  load  position 
overshoot  was  of  interest  is  given  in  Appendix  0.  Approximately  1400 
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phase  trajectories  were  solved  using  this  mode  of  operation.  Each 
solution  required  an  average  of  three  minutes  of  computer  time. 

For  the  purpose  of  general  examination  of  the  problem  it  was 
assumed  that  P  =  1,  =  1,  0\\=  f  -^L  =  1-  Specific  solutions 

were  made  with  p  ^  1  (  1  to  determine  scaling  effects.  The  re¬ 

sults  of  these  solutions  are  presented  in  Section  6,  Application  of 
Results. 

The  program  can  also  be  made  to  solve  linear  systems  by  setting 

A  =  0. 


The  major  variable  parameters  used  in  this  investigation  were: 

^  =  o.i,  0.2,  0.3,  0.4,  0.5,  0.6,  o.8,  l.o 

o  ,  0.2,  0.4,  0.6,  0.8,  l.o 

"^T 

9_5-=  0.1,  0.2,  0.5,  Q.8,  0.9 

oTU,  0.5,  0.2,  0.1 

A  =0.3  radians 


e  =  0,  0.6,  0.8,  1.0 
Additional  parameters  of 


0.04,  0.1,  0.9;  and  A  =  0.01,  0.03,  0.01,  0.15  radians 
were  used  in  certain  instances  to  examine  particular  characteristics  of 
the  response. 

The  program  was  limited  somewhat  in  that  the  value  of  P^v-  =  qq 
and  zero  were  exluded  due  to  generation  of  undefined  mathematical 
quantities. 

The  value  of  backlash  was  made  abnormally  large,  A  =0.3  radians, 
to  allow  an  easier  interpretation  of  the  non-linear  response.  Since 
the  influence  of  A  was  linear,  this  caused  no  inaccuracies.  All 
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graphs  for  limit  cycle  size  were  plotted  with  /\  =  0.3  radians. 

The  effect  of  values  for  /\  other  than  0.3  is  also  shown  in  graphi¬ 
cal  form. 
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5-  Results  and  Discussion 


Two  types  of  system  response  were  observed.  The  first  type  was 
a  convergence  to  zero  error,  which  was  a  case  of  static  and  dynamic 
equilibrium  of  the  system.  This  is  a  characteristic  of  a  linear 
second  order  servomechanism.  The  second  type  of  system  response 
observed  was  the  divergence  or  convergence  to  a  state  of  dynamic  equi¬ 
librium.  This  second  type  of  response  is  termed  a  limit  cycle  and  is 
characterized  by  cyclic  travel  of  the  load  through  the  same  points 
of  the  phase  trajectory.  When  the  second  mode  of  computer  operation 
was  used,  both  the  maximum  load  position  per  cycle  and  the  period 
of  the  recorded  position  were  analyzed  to  determine  if  the  character¬ 
istic  response  was  a  limit  cycle.  Exact  repetition  of  cyclic  values 
was  impossible  due  to  the  numerical  methods  used  and  the  computer 
round-off  error  for  the  output  routine.  Typical  examples  of  the 
limit  cycle  and  no  limit  cycle  response  are  given  in  Appendices 
D  1  and  D  2. 

No  case  of  dynamic  instability  or  divergence  of  the  system  with¬ 
out  bound  was  observed  in  the  investigation.  Intuition  might  lead 
to  this  same  conclusion  when  the  following  unique  characteristics  of 
this  system  are  pointed  out:  (a)  As  a  maximum  limit,  energy  was 
conserved  at  the  instant  of  impact,  (b)  Energy  was  supplied  to  the 
system  only  in  finite  amounts  and  at  a  finite  rate  by  the  error 
detector.  (c)  Energy  sinks  were  present  in  the  system  frictional 
elements  while  energy  storage  units  were  present  in  both  the  motor 
and  load.  (d)  No  undefined  limits  appeared  in  the  applied  physical 
equations, 


O  vV\ 


0  and  excluded. 
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The  results  of  the  investigation  are  offered  in  the  form  of 

charts  for  parameter  areas  of  limit  cycle  existence  and  nonexistence 

in  Fig.  1  for  e  =  0.,  0.6,  0.8  and  Fig.  2  for  e  =  1.0.  The  para- 
Z  r  2. 

meters  P  ti_  and  p  3l  are  designated  abscissa  and  ordinate 

respectively.  The  parameter  points  of  examination  for  which  solutions 
were  obtained  are  circled.  The  zone  enclosed  by  a  particular  ^  line 
is  largely  an  interpretation  by  the  authors  of  the  results  of  the  in¬ 
vestigation,  and  shows  the  area  in  which  no  limit  cycle  existed.  Due 
to  the  large  change  of  the  variable  parameter  values  between  solutions, 


the  charts  are  not  exact.  All  circled  intersections  interior  to  the 


s 


line  indicate  parameters  which  resulted  in  no  limit  cycles,  e.g.. 
Fig.  1  for  e  =  0,  0.6,  0.8  at  a  ^  =  0.4,  the  values  of  = 


Ft 


0.6,  0.8  at  a  ^  =9*0  and  the  value  of  ^  it  =  0.6  at  a 

0  wt 

=  4.0  resulted  in  no  limit  cycle.  Similarly,  for  =  0.5, 


3  Wt 


the  values  of  =  0.4,  0.6,  and  a  O  =1.0  as  well  as  all 

TV 

other  parameter  intersections  interior  to  the  =0.5  curves,  re¬ 

sulted  in  no  limit  cycles. 

All  circled  parameter  intersections  exterior  to  particular  ^  lines 
define  paramters  which  were  examined  and  resulted  in  a  steady  state 
limit  cycle.  The  area  to  the  right  of  the  single  ^  =1.0  line 
designates  the  parameters  for  which  no  limit  cycles  were  observed. 

Limit  cycles  were  observed  for  all  ^ ^  0.5  for  the  parameters  examined. 

The  investigation  was  made  with  ^  =0.5  radians;  however,  re¬ 

sults  of  the  investigation  of  variation  of  limit  cycle  size  indicate 
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that  existence  of  the  limit  cycle  is  independent  of  the  size  of  back¬ 
lash  in  a  nonlinear  system  for  a  A.  greater  than  zero.  Variation  of 
A  had  the  effect  of  a  linear  variation  of  limit  cycle  size  only. 

The  zones  depicted  on  the  existence  charts  are  approximately 
represented  by  straight  lines;  however  there  is  reason  to  believe  that 
a  complete  examination  would  show  a  curvature  in  the  lines.  This 
observation  will  be  elaborated  upon  later  when  the  relationship  between 
the  existence  charts  and  the  figures  of  limit  cycle  size  is  discussed. 

From  Fig.  1  it  is  seen  that  there  is  only  a  very  small  difference 


in  the  limit  cycle  existence  zones  for  plastic  impact  (e  =  0)  and  the 
intermediate  values  of  the  coefficient  of  restitution,  e  =  0.6  and  0.8. 
When  comparing  the  elastic  impact  (e  =  l)  in  Fig.  2  to  the  plastic  and 
intermediate  cases  in  Fig.  1,  the  most  obvious  point  of  comparison  is 
the  drastic  reduction  of  the  areas  in  which  limit  cycles  did  not  result 
for  e  =  1.0.  However,  it  may  be  generally  observed  that  with  any  e 
held  constant,  the  area  in  which  a  limit  cycle  could  not  be  obtained 
increased  with  the  increase  of  system  ^  . 

For  all  values  of  e,  there  is  an  apparent  common  enclosed  zone 

centroid  which  is  described  by  the  system  parameters  yy  f ^yF^  =0.5 

Vi¬ 
and  p  ^L/^m  =  1*0'  These  parameters  result  in  the  system  time  con¬ 
stants  =  since  £0n  =  1  ani^  =  1  were  chosen  for  the 

investigation.  There  is  also  an  apparent  symmetry  about  the  centroid 
of  opposite  quadrants  which  can  be  related  to  the  fact  that  since 
(j)=  1  and  p  =  1,  the  points  diametrically  opposite  across  the 


centroid  of  the  figure  (using  the  decimal  values  of  the  parameter  inter¬ 
section)  represent  an  exact  exchange  of  component  time  constants.  (See 
Fig.  2  for  example).  Two  diametrically  opposite  points  may  have 
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exactly  the  same  motor  and  load  time  constants;  however,  they  are 
descriptive  of  two  different  systems  since  the  friction  and  inertia 
distributions  are  different  for  the  two  points. 

Figs.  5  through  26  present  the  maximum  positive  error  in  radians 
of  the  limit  cycle  for  a  unit  step  input  and  with  backlash  equal  to 
0.3  radians.  The  use  of  the  charts  is  explained  in  Section  6,  Appli- 
cation  of  Results.  Each  figure  is  presented  with  p  f^yF  and  e 


constant.  The  independent  parameter 


,J  is  indicated  as  abscissa 


?  V- 

and  curves  of  constant  3  are  plotted  with  limit  cycle  error  as  the 

-t 

ordinate.  The  results  were  presented  on  a  log  log  plot  since  the 
observed  limit  cycle  errors  encompassed  a  range  from  0.01  to  1.0 
radians,  and  the  appearance  and  disappearance  of  limit  cycles  and 
the  change  of  limit  cycle  error  was  more  readily  apparent  on  the  log 
log  form. 

In  general  as  e  increases,  the  size  of  limit  cycle  decreases,  all 

other  parameters  remaining  constant.  However  this  is  not  without  a 

few  exceptions; the  generalization  is  most  accurate  for  large  size  limit 

cycles.  It  may  be  noted  that  when  viewing  the  individual  charts,  if  a 

k 

limit  cycle  exists  at  very  low  P  J  /j  ratio  there  is  a  tendency  for 

\  L  M 

the  same  system  to  have  a  smaller  limit  cycle  for  higher  P  J  /j  ratio. 

\  L  m 

In  addition,  for  these  same  conditions,  if  limit  cycles  exist  for 

several  values  of  ^  ,  there  is  a  convergence  toward  the  same  size  of 

limit  cycles  for  each  ^  as  pj/j  increases.  At  the  value 

\  L  m 

"i-  . 

p  f^/F^  =  0,  a  limit  cycle  will  always  occur  for  any  value  of  e  or 

Viewing  the  charts  in  sequence,  with  either  e  constant  or  <p  f^/F 
constant,  it  is  noted  that  there  is  continuity  of  pattern  flow  between 
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adjacent  charts.  This  continuity  between  charts  appears  to  form  a 

symmetrical  pattern  of  constant  lines  which  can  be  interpolated  for 

L 

intermediate  values.  When  the  figures  for  p  f  /F  =  0.4  and  0.6 

V  L'  t 

are  analyzed,  this  symmetrical  pattern  is  observed  in  each  chart  about 

p  J  /J  a  l.  The  limit  cycle  decreases  and  disappears  at  the 
V  L  m  t 

The  charts  of  f)f/F  =  0.4  and  0.6 
V  1/  T 


opposite  extremes  of  p  J  /j 

1  L  m 

are  symmetrical  to  each  other,  while  ^j/'^T  =  0*2  and  0.8  are 
similar.  At  the  extreme  values  of  ^  f  /F^,  =  0  and  1.0,  no  symmetrical 
comparison  can  be  made.  An  additional  point  of  symmetry  about 
^  f^/F^  ~  can  be  visualized  for  all  values  of  e.  These  para¬ 

meters  for  symmetry  are  noted  in  the  discussion  of  the  existence  charts . 
The  irregularities  in  the  curves  of  medium  values  of  ^  can  not  be 
explained  by  the  authors;  however,  it  may  be  mentioned  again  that  the 
points  of  symmetry  are  characterized  by  exact  exchange  of  component 
time  constants. 

The  noted  symmetry  of  the  limit  cycle  size  charts  shows  a  direct 
relationship  to  the  symmetry  found  in  the  existence  charts.  Since  the 
existence  charts  show  contours  of  constant  ^  for  a  surface  of  zero 
limit  cycle  size,  one  is  led  to  the  possibility  of  a  three  dimensional 
presentation  to  show  a  variation  of  limit  cycle  size  for  the  same  para¬ 
meters  used  in  the  existence  charts.  Such  a  three  dimensional  figure 
for  constant  e  using  ^  f^/F^,  and  ^  as  the  axes  is  shown 

in  Fig.  27 ♦  Other  surfaces  could  be  added  above  this  zero  limit  cycle 
"floor"  to  show  the  variation  of  a  given  limit  cycle  size  by  contours 
of  ^  .  The  ultimate  result  of  this  line  of  investigation  would  be 
the  evaluation  of  the  poles  and  zeros  for  a  confoimal  mapping  presen¬ 
tation.  Unfortunately,  time  did  not  permit  the  authors  to  investigate 
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this  avenue  flirt  her. 

The  results  of  the  investigation  of  the  dependence  of  limit  cycle 
size  on  the  amount  of  backlash  are  presented  in  Figs.  28  and  29  using 
coefficients  of  restitution  of  0  and  1.0  at  a  ^  =0.6.  A  thorough 
investigation  was  made  at  this  ^  since  it  is  an  average  system 
parameter  often  used  in  design  work.  The  limit  cycle  size  was  found 
to  vary  directly  with  backlash  size  without  exception  between  the 
values  of  0.3  and  0.01  radians.  When  limit  qycles  existed  at  ^  =  0.3, 
there  was  no  disappearance  of  limit  cycle  with  decreasing  magnitude  of 
A.  .  Where  no  limit  cycle  existed,  at  /\  =  0.3,  limit  cycles  con¬ 
tinued  to  be  non-existent  throughout  the  range  of  examined. 

Several  solutions  of  problems  with  =  0.6  and  e  =  0.6  &0.8  and  = 
0.1&e=1.0&0  with  various  values,  proved  the  linear  variation 
of  limit  cycle  size  with  backlash  size  for  all  values  of  e. 

When  a  servomechanism  under  steady  state  conditions  is  disturbed 
by  a  small  perturbation,  it  will  return  to  the  original  steady  state 
conditions.  However,  the  transient  response  due  to  the  perturbation 
may  be  quite  different  from  the  transient  response  which  occured  due  to 
the  original  signal  input.  The  concepts  of  energy  and  static  and  dy¬ 
namic  equilibrium  mentioned  in  Section  1  may  be  better  understood  as  a 
result  of  the  perturbation  technique  if  the  usual  initial  conditions  of 
4  (0)  =  4  (0)  =  0  are  considered  as  a  disturbance  from  the  stability 
conditions  predetermined  by  selection  of  system  parameters  and  signal 
input. 

This  same  small  perturbation  technique  was  simulated  mathematically 
using  the  digital  computer  to  obtain  several  solutions.  The  cases 
examined  and  their  initial  conditions  were: 
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a)  ^  =  0.6,  ^JL/Jm  =  0.8/0.2,  p"fL/FT  =  0,  e  =  0,  A=0.3 

-6cc  (0)  =  0  ;  -6  c  (0)  =  1.2 

4  m  (0)  =  0  ;  0m  (0)  =  1.0 


system  separated. 


system  combined 


c)  '§=0.4,  P'j/J  =  0.5/0. 5, 

\  L  m 

c  (0)  =  0  j  &  c  (0)  =  1.005 

* 

■0-  m  (0)  =  0  ^-m  =  1.00 

•J 


system  separated 

In  the  three  cases  examined,  the  same  limit  cycle  was  obtained, 
both  in  size  and  period,  as  that  resulting  from  the  normally  used 


between  the  perturbation  response  and  the  original  limit  cycle  size 


for  the  three  cases  was  18  x  10-^  radians.  Two  types  of  transient 


response  were  observed.  These  were  convergence  to  limit  cycle  in  the 
first  two  cases  and  divergence  to  limit  cycle  in  the  last  case. 

The  results  obtained  by  the  small  perturbation  technique  supported 
the  validity  of  the  analysis  of  the  problem  by  digital  methods.  This 
result  also  concurs  with  the  findings  for  the  case  of  e  =  0  by  Knoll 
and  Narud,  Ref.  d,  in  that  the  steady  state  limit  cycle  size  is  completely 
specified  by  the  choice  of  system  parameters  and  independent  of  system 
initial  conditions. 

Mention  is  made  here  that  the  computer  program  could  not  be  run 
normally  for  the  perturbation  method.  The  computer  was  stopped  after 
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the  print  out  of  system  constants,  then  the  initial  conditions  of  the 
servo  system  were  inserted  manually.  The  program  was  then  restarted 

at  the  desired  point  for  computing,  either  as  a  combined  system  or 
with  the  motor  and  load  separated.  The  program  was  then  allowed  to 
operate  to  the  normal  completion  of  the  solution. 

The  work  of  Ref.  d  reported  that  the  system  transient  response 
for  the  case  e  =  0  arrived  at  a  steady  state  limit  cycle  by  two 
possible  phase  trajectories.  One  of  these  phase  trajectories  was 
smooth  convergence  to  the  limit  cycle,  while  the  second  was  a  converging 
overshoot  toward  the  inside  of  the  limit  cycle  followed  by  a  divergence 
to  the  limit  cycle.  From  the  results  of  this  investigation  for  the 
cases  e  =  0,  0.6,  0.8  and  1.0,  similar  transient  indications  were  ob¬ 
tained. 

Several  computer  solutions  were  accomplished  for  the  cases  =0.3 
radians,  e  =  0,  1.0,  p  =  1.0,  n  =  1.0,  K  =  J  f  1.0.  In  addition 
various  combinations  of  p  /  1.0  and  ^  n  /  1.0  were  obtained.  The 
purpose  of  these  solutions  were: 

1.  To  determine  the  general  applicability  of  the  computer  program. 

2.  To  verify  the  scaling  developed  in  Section  2,  Development  of  Equations 

3.  To  determine  the  applicability  of  the  results  of  this  thesis. 

The  following  results  were  determined: 

A.  Two  solutions  were  obtained  (e  =  0  &  1.0)  for  the  case  ^  =  0.4, 


P  =  1.0  U)  n  =  1.0,  J  =  2.0,  K  =  2.0,  W  L=-  =  0.4,  4.0. 

The  results  of  these  solutions  were  the  same  (time,  load  velocity,  load 

position,  and  limit  cycle  size)  as  the  solutions  obtained  for  Vj  <-  =0.4 
■? 

D  3^  =  4.0,  ^  =  0.4  using  the  nondimensional  parameters  (j)  n  =1.0, 

'  3T  •'vi 

p  =  1.0,  J  =  1.0,  K  =  1.0. 
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B.  Solutions  were  obtained  for  e  =  0  &  1.0  for  the  case  =  0.4, 

=  1.0,  00  n  =  2.0,  J  =  1.0,  K  =  4.0,  =  0.4,  J?Ti_  =  4.0. 

Ft  '  Twi 

The  results  of  these  solutions  yielded  the  same  size  maximum  overshoots 


e 


§  =  0.4, 


and  limit  cycle  size  as  the  nondimensional  solutions  for 

Pjh_=  0.4,  ^  .0 1  -  4.0.  The  times  of  maximum  overshoots  in 

'  Ft  3"  vv\ 

this  case  were  one~half  the  times  of  occurence  of  the  same  overshoots 


for  the  nondimensional  solution.  The  velocity  at  a  maximum  overshoot 


was  twice  the  value  of  velocity  obtained  for  the  same  overshoot  in  the 
nondimensional  solution.  This  same  comparison  can  be  made  for  Cases  A 
and  B.  The  results  of  these  bnd  other  solutions  prove  the  validity 


formulae  used  in  Section  2,  Development  of  Equations,  concerning  the 
transformation  to  the  coordinate  system.  (*) 

C.  Solutions  were  obtained  for  the  case  e  =  0  &  1.0,  'S  =  0.4,  p  =  0.5, 

£l)n  =  1.0,  4j_  =  0.8,  Tl_  =  4.0,  =  0.2,  =1.0.  The 

f~T  T  Ft  t  wv 

results  of  these  solutions  were  the  same  (time,  load  velocity,  maximum 


overshoots  and  limit  cycle  size)  as  for  the  nondimensional  solutions 

=  o.4,  P  =  1.0,  6d n  =  1.0,  =  0.2,  P1_7u=  1.0. 

v  Ft  ' 

D.  Solutions  were  obtained  for  the  case  e  =  0  &  1.0,  ^  =  0.4, 

P  =  0.5,  lj)  n  =  2.0,  =  0.8,  =  4.0,  (£lt-  =0.2,  ^  Jl  =  1.0. 

Pf  TTyv'x  FT  "3  yv\ 

The  results  of  these  solutions  yielded  the  same  size  maximum  overshoots 


and  limit  cycle  size  as  the  nondimensional  solutions  for  =  0.4, 


& 


)^L=  0.2,  P  =  1.0.  The  times  and  velocities  in  this  case  were 

Ft  ,  .  Sw 
scaled  by  60  n  as  in  Case  B. 


Several  phase  trajectories  using  e  =  0  were  compared  with  those 


obtained  by  New,  Ref.  e.  Generally  good  correlation  was  obtained  on 
the  value  of  limit  cycle  error  and  in  the  transient  response.  Small 


disagreements  were  expected,  resulting  primarily  from  differences  in  the 
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programming  procedures  and  analysis  of  results  by  read-out  methods. 


Sample  phase  trajectories  for  the  cases  of  e  =  0  &  1.0  are  pre¬ 
sented  in  Figs.  30  and  31.  Although  the  transient  analysis  of  the 
system  is  beyond  the  scope  of  this  work,  these  figures  may  aid  in 
understanding  the  physical  problem. 

Of  academic  interest  was  the  accidental  operation  of  the  computer 


program  for  several  solutions  using 


Several  responses  were 


observed.  One  was  a  nearly  pure  oscillatory  system  from  the  first 
overshoot  and  thereafter  another  was  a  very  slow  convergence  or  slow 
divergence  of  the  oscillations  after  the  first  overshoot.  The  con¬ 
vergence  or  divergence  probably  resulted  from  errors  generated  in 
numerical  iterations  and  approximations. 
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6.  Application  of  Results 

In  this  section,  the  results  of  this  thesis  will  be  applied  to 
hypothetical  physical  systems.  Methods  will  be  indicated  by  which  physical 
systems  may  be  analyzed  or  designed.  Attention  is  called  to  the  assump= 
tions  stated  in  Section  2  under  which  this  study  was  made. 

It  is  desirable  to  find  an  equation  relating  certain  system 
variables,  therefore  the  equations: 


a) 


^ 

vi  yyy  -V  T  L 

a  v\  d 


-  P 


k 


(2) 


tp  f 

P; 


(  m  v  n 
P^L  J 


P 


will  be  manipulated  to  the  form 

<5)  ^  =  -1  fit  - 

0  w\ 

or 


+  P  3l 


DJl 

T*1 


M 


■h  \ 


4-  ^ 

1  w\| 


VTtn  \|  \  +  p"2t 

vHT  \nr  '  dm 


C  -  J-  £ ,  \  f  "-1  4-1 

1  -1 

P 

00 

^  x  LpMu  _ 

-  3  vn  C  4-  ^  vW  ^ 

'4 


In  view  of  the  form  of  Equation  (4)  for  ^  and  recalling  that  the 
variable  parameters  used  on  the  figures  were  P  3-b  and  P  ^i=.  = 

^  r  \  '  "Ft 

V  •“  l.  =  r  ..  ,  \  it  is  seen  that  the  expression  for  -P 

r - ^  1  O 

T  \  P  V  t_  p1-  r  1- 

in  fact  includes  the  other  two  variable  parameters. 
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Fran  discussion  concerning  the  effects  system  natural  frequency, 
it  will  be  recalled  that  the  value  of  overshoot  and  hence  limit  cycle 
are  not  affected  by  CO  n;  however,  the  time  and  velocity  are  scaled 
for  n  ^  1.0.  Thus  the  transient  response  characteristics  may 

be  related  to  the  same  parameters  used  in  these  applications. 

If  a  motor  and  load  are  selected  for  a  system  and  the  limitation 
is  imposed  that  system  operation  shall  not  result  in  a  limit  cycle, 
Equation  (4)  and  the  limit  cycle  existence  charts  may  be  used  to 
determine  the  unspecified  system  parameters  ,  K  and  ^  .  The 

coefficient  of  restitution  for  the  proposed  gear  train  material  may 
be  determined  by  the  use  of  the  equation  in  Section  2. 

To  proceed,  select  a  value  of  K  and  (p  ,  then  determine  the  value 
of  ^  from  Equation  (4).  Enter  the  limit  cycle  existence  chart  for 

the  value  of  e,  system  and  the  parameters  P  ^  *-  and  — — ^ -  . 

-*-1 


pz(i 


Determine  if  the  selected  parameters  describe  a  point  of  no  limit 
cycle.  Since  limit  cycle  existence  is  independent  of  backlash  ( f  0 
excluded),  backlash  is  not  a  variable  parameter.  Successive  trial 
values  of  the  unspecified  parameters  may  be  required  to  meet  the 
limitation  of  no  limit  cycle.  If  it  is  also  required  that  certain 
transient  response  characteristics  be  met,  these  conditions  could  be 
examined  by  a  similar  method. 

In  the  event  that  no  acceptable  parameters  are  found  to  satisfy 
the  system  requirements  for  no  limit  cycle,  the  figures  for  limit  cycle 
size  may  be  examined  to  obtain  a  minimum  size  limit  cycle.  The  method 
of  solution  is  similar  to  that  previously  described. 

Example:  The  material  intended  for  use  in  the  gears  has  an  e  =  0.6. 
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The  motor  and  load  have  the  parameters  fin  =  0.64,  Jm  =  0.25,  fT  =  0.64, 

.L 

=  1.0  and  it  desired  that  the  system  operate  with  no  limit  cycle. 

If  the  values  are  selected  K  =  4.0,  ^  =  0-5,  the  values  =  1.0 

4m  =4.0.  - L  ~  ~  -■  ^  ™ 


and  i  M 


e'f 


=  0.2  result. 


.  5 


f  ^  +■  1 

Solve  Equation  (4)  for  t  -  .5(.64)(5+l)  _  _ 

=  0.4.  Enter  the  chart  for  limit  cycle  existence.  Fig.  1, 
e  »  0.6,  'C  -  0.4,  — - J -  -  0.2  .  C?  31  •  1.0 

*  Sin-'  '  ~ 


,  thus 

2. 


p*C< 


It  can  be  seen  that  limi-t  cycle  will  exist  for  these  parameters. 
Since  K  is  most  easily  varied  in  equation  (4),  a  new  K  is  selected 
at  1.0  in  an  effort  to  find  parameters  for  no  limit  cycle.  It  is  seen 


that 


t 


K  =  U 


t 


K--  H.o 


=  0,4  (?) 


with  all  other  parameters  constant,  thus  ^  =0.8  Enter  Fig.  1 

rtth  e=0.6  <).*  ^  ,J  ()  J 


0  wi 


J 


-  0,1 


fi 


_J21  +  l 

These  parameters  describe  system  operation  with  no  limit  cycle  and  are 
well  within  the  area  limits.  The  transient  response  of  this  system 
may  be  undesirable  and  various  combinations  of  p  and  K  would  have  to 
be  tested  to  satisfy  the  additional  requirements  of  the  problem.  It 
will  be  restated  that  the  existence  areas  of  Figs.  1  and  2  are  not  exact. 

The  figures  of  limit  cycle  size  may  be  used  in  the  same  manner  as 
the  existence  charts.  The  same  parameters  of  the  first  case  are 
assumed:  K  =  4.0,  p  =  0.5>  ^  =  1.0, 

O  VY\ 

54 


ir  +i 

ru 


=  0.2 


K  -  0.4 


e  =  0.6  A  =0.1  radians.  Enter  Fig.  10  with  the  parameters  given. 

The  ordinate,  Limit  Cycle  Error  (  0. 3  \  =  0.13  radians,  is  obtained. 

V  A  / 

The  maximum  positive  error  in  the  limit  cycle  =  0.13  (O.l)  =  0.0^3^ 

radians.  Enter  the  same  figure  for  e  =  0.6,  ^  =0.8  )  (  Ju  =  1.0 

I  .  “ 

i  »  =  0.2  ,  =  0.5  radians;  this  is  the  second  case  of  exami- 

pFu  1  ; 

nation  for  existence.  No  \T  =  0.8  curve  exists;  hence  no  limit  cycle 
exists  for  this  choice  of  parameters  for  any  nominal  value  of  backlash. 


) 
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7.  Conclusions 


Within  the  limitations  of  the  assumptions  made  and  the  scope  of 
this  thesis,  it  is  concluded  that: 

A.  Two  types  of  system  steady  state  response  result.  The  type  of  re¬ 
sponse  is  independent  of  backlash  size,  system  natural  frequency  and 
initial  conditions.  The  types  of  response  are  convergence  to  signaled 
input  position  and  convergence  to  a  limit  cycle. 

B.  If  a  limit  cycle  exists  in  the  steady  state  response  of  a  system, 
the  size  of  the  limit  cycle  varies  directly  with  the  size  of  backlash, 
all  other  parameters  constant*. 

C.  The  greatest  change  in  system  steady  state  response  occurred  for 
the  case  when  system  rotational  energy  is  conserved  (e  =  1.0). 

D.  A  point  of  symmetry  for  the  type  of  system  steady  state  response 
exists  with  equal  motor  and  load  time  constants.  For  medium  to  high 
values  of  system  damping  coefficient,  symmetry  with  respect  to  limit 
cycle  size  exists  for  an  exchange  of  motor  and  load  time  constants. 

E.  The  results  of  this  thesis  are  presented  in  a  form  which  may  aid 
in  system  design  or  analysis. 


36 


8.  Recommendations 

From  the  results  of  this  thesis  it  can  be  seen  that  there  remain 
many  fertile  areas  for  further  investigation.  For  studies  of  this 
type  it  is  strongly  recommended  that  the  digital  computer  be  used 
because  of  its  speed  and  versatility.  The  limits  of  investigation 
depend  to  a  great  extent  on  the  ingenuity  of  the  programmer.  Some 
recommended,  areas  of  future  investigation  are: 

a)  More  precise  definition  of  the  limit  cycle  existence  zones 

and  a  more  thorough  coverage  of  the  limit  cycle  size  of  similar  systems 
having  the  same  type  of  parameters.  A  three-dimensional  presentation 
as  mentioned  in  the  discussion  could  be  obtained. 

b)  Include  torsion  and  other  deformations  of  the  mechanical  system 
in  the  physical  equations  of  the  system  to  determine  their  affect  on 
the  system  response. 

c)  An  analysis  of  the  limit  cycle  existence  zone  and  transient 
response  of  a  system  to  determine  if  there  is  a  relationship  between 
the  two.  An  analytical  expression  for  the  possibility  of  limit  cycle 
existence  might  be  developed  in  terms  of  system  parameters. 

d)  Examine  the  transient  and  frequency  response  of  the  system  due 

to  various  input  functions.  This  would  require  minor  modifications 

to  the  computation  section  of  the  existing  computer  program. 

e)  Since  the  physically  realizable  quantities  of  position,  velocity 
and  acceleration  of  both  the  motor  and  load  are  available  as  computed 
quantities,  the  study  of  the  compensation  problem  and  its  optimization 
are  obvious  areas  for  further  consideration. 

f)  The  general  study  of  other  higher  order  nonlinear  systems  is 
suitable  for  digital  computer  programming  since  most  numerical  inte- 
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gration  methods  will  handle  systems  having  equations  of  any  order. 
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APPENDIX  A 

COMP  U  TA  T/OSV  FLOW  £>/A  G-PA  M 

for 

CDC  /604  D/G-/TAL  COMPUTER 


7Z 


from  MOTOR.  3Q 
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from  S9,  GO 


from  S£,6 


from  4*3 


83 


EXIT’  _ 

84 


ref  urn  to  progrmtr 


74 


APPENDIX  A  (continued) 


1. 


2. 


3- 

k. 

5- 

6 . 


7- 

8. 

9- 


10. 


ll. 


12. 


13. 

14. 


EXPLANATION  OF  FLOW  DIAGRAM  BLOCK  SYMBOLS 
Clear  fault  stop. 

Print  last  overshoot  stored  in  LIMIT3UF  (symbol  888). 

Print  out  position  in  phase  trajectory  of  fault  stop 

(symbol  666). 

End  of  file. 

Actuate  fault  stop. 

Load  zeta  from  Index  1  (b  ). 

1 

Print  last  overshoot  stored  on  LIMITBUF  (symbol  888) . 

Print  position  in  phase  trajectory  where  it  was  manually 
stopped  (symbol  1961). 

Stop  if  Zeta  -  0. 

2JO).  -  a 

Is  zeta  -  0,1? 

Is  zeta  less  than  0.1? 

Set  up  computation  time  of  0.01  seconds  for  problem  and 
set  up  for  print  at  every  0.1  seconds  when  zeta  is  greater 
than  0.1. 

Set  up  computation  time  of  0.00^  seconds  for  problem  and.  set 
up  for  print  at  every  0.1  seconds  when  zeta  is  equal  to  or 
less  than  0.1. 

Set  up  limit  cycle  print-out  routine  for  20  overshoots  and 
average  the  last  8. 

Set  up  number  of  limit  cycle  print-outs  and  averaging  routine 
depending  upon  the  size  of  zeta: 
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f 

Overshoots 

Start 

Averaging 

Divide 
by  * 

Comp. 

Time 

^  0.8 

6 

4 

2.0 

0.01 

0.6 

7 

5 

2.0 

0.01 

0.5 

8 

5 

3.0 

0.01 

0.4 

10 

7 

3-0 

0.01 

0.3 

14 

10 

4.0 

0.01 

0.2 

20 

12 

O 

CO 

0.01 

<  0.1 

17 

12 

5.0 

0.004 

*The  last  row  printed  has  this  number  printed  on  the  far  left 
if  the  problem  has  a  normal  solution.  This  number  indicates 
how  many  overshoots  were  averaged.  The  average  value  has  1.0 
subtracted  from  it  and  is  printed  in  the  far  right  of  the  same 
last  row. 

15 .  Page  re  j  ect . 

16.  Clear  buffers  and  set  up  for  new  problem. 

17*  Restitution,  e,  from  Index  3  (b^,). 

Backlash,  ,  from  Index  4  (b^). 

J  rr,  from  Index  5  (b,_ )  • 

5 

Jt  “  Jw  K  -  Jtn 

jjt-  from  Index  6  (bg). 

(yfTy -  f*-  j  f!yfL  -  & 
v°2fL  =  &  j  -  0 

A  r 

18.  Have  all  values  been  cycled  through? 

19.  Set  STOPNOW  after  this  run  to  stop  absolutely. 


=  0? 


20.  Backlash,  ^  , 


21.  Print  consecutive  number  of  this  parameter  run  and  all  the 
parameters  used: 

£j°6J»=A  J'Yjl  Vft  & 

JL  JmV9*!  frn  fL  /°2 

k  K/j^c  Vj.-o  / 

22.  Set  program  so  always  acts  as  combined  system. 

23.  Print  consecutive  number  of  this  parameter  run  and  selected 
parameters  for  linear  system: 

Zf6>*=A  K  GO*? 


/°  /° 


2k.  Solve  the  second  order  equation  by  Runge-Kutta-Gill  numerical 

integration  every  0.01  seconds  problem  time  when  zeta  is  greater 

than  0.1,  otherwise  every  0.004  seconds.  Initial  conditions 

•  •« 

are  1«0,  ©c  =  ~  ®  an<^  =  °r  Pre“ 

viously  computed  point.  Next  point  is  computed  by  four 
iterations  of: 

0C~  ^ h  (©r  ~  ©c)  “  &  jf*  CO  h  0c 

25.  Set  EXIT  from  print  routine. 

26.  Load  time,  ac  and  in  print  buffer. 

27.  Go  to  LDdTCYC  to  print  only  the  overshoot.  Otherwise  print 
the  phase  trajectory  points  for  each  0.1  seconds  up  to  a  maxi¬ 
mum  of  12  cycles. 

28.  Store  0  for  £)  and  ©  m 
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29. 


30. 


31. 


32. 


33. 


34. 


35- 

36. 


37- 


38. 


39. 

40. 

41. 

42. 

43. 

44. 
45- 
46. 


Is  problem  time  0.1  sec.? 

Is  problem  time  greater  than  0.1  sec.? 

o 

Print  the  number  of  point  computed,  time,  ©c  and  ©c  . 
EXIT  (automatically  set  for  desired  jump-out.) 

O®  ° 


Store  time,  ©C  ’  0c  and  ©c  for  motor  and  load  initial 

conditions. 

Is  slope  Nt  =  N_? 

-L'  S 

(©R-G>c)d0n_3  __  * 

- sr-1 —  eLy  ^>v\T-  J-h-  =  positive  value? 

J  L 

Is  Q  negative? 

Store  ©c-  A  irr 


/° 


Solve 


©c  = 


for  ©c  and  ec 


by  Runge-Kutta-Gill. 


for  and  0*  by  Runge-Kutta-Gill 

0L  ”  y0  ^  o  ? 

0  L  ~ /°  Q  A  =  negative? 


Set  EXIT  from  print  routine. 

Set  EXIT  from  print  routine. 

Load  time,  and  Q  c  , 

o 

Load  and  - 

Print  number  of  point  computed  and  stored  values. 
EXIT  (automatically  set  for  desired  jump-out). 
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4-7  •  Solve 

©  / 

Oc  = 


I  KKV 


<9m  (i  +  e)  -t  ©c  (/?  -A  -  e) 


©m 


©*  s 


&  -  e(/<?erh  +■  ©c) 


o  / 

-h 


7-  (©*  -  ©0 


4-8.  Set  EXIT  from  print  routine. 

© 

49.  Is  ©c  positive? 

© 

50.  Set  up  to  count  one  phase  trajectory  when  ©c  is  negative. 

51.  Is  0C  _>  1.0? 

52.  Have  12  cycles  of  phase  trajectory  "been  completed? 

53*  Arrange  for  no  counting  until  &c  is  positive  again. 

*  A 

5^.  Is  ©c  negative? 

55 •  Set  up  —  (jr  ®  n*  “  S>c)  . 

*  © 

56.  Is  "*  greater  than  preselected  £  ? 

57*  Set  up  combined  system. 

58.  Is  ©c  positive? 

© 

59.  Is  ©c  negative? 

60.  Is  G>c-  l  positive? 

61.  Is  —  I  negative? 

62.  Is  0c“  I  negative? 

63-  ©c-^@m=0? 

64.  Oc-fBn=0l 

65.  Is  <3^  “  ^  positive? 

66.  Is  @c  ~yO  Prepositive? 
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67. 


68. 

69- 

70. 

71. 

72. 


73. 

74. 


75- 

76. 


77- 


78. 


79- 

80. 

81. 

82. 

83. 

84. 

85. 

86. 

87. 


Is  Slope  =  Kg 

C I  -  Sc)  U)n 


O 

©C 


-  + 


fl 


Modify  FLOATEST  to  keep  system  combined. 


when  ©e  is  negative. 

Is  0^  —  Syy\  positive? 

Store  ©L.  ~  A  in  ^ 

Store  0C  in  0  ^ 


=  positive  value? 

Takes  off  modification 


Is  ©C  ^  I  negative? 

0C  =  oi 

Sc  ~  fOv*=  0? 

Is  ©e  “ ^0  ©^positive? 

Is  0C  -y/9  ©**  positive? 

Modify  FLOATEST  to  keep  system  combined.  Takes  off  modification 

c 

when  Sc  positive. 

Have  ll/2  minutes  of  real  time  elapsed  since  the  last  print-out? 
Is  @c  <  1.0? 

Set  up  to  print  when  in  the  fourth  quadrant. 

Check  if  in  the  fourth  quadrant. 

Is  the  new  ©e  greater  than  the  previous  Oc  ? 

* 

Store  the  new  time,  ©c  and  ©<:• 

EXIT  (automatically  set  for  desired  jump-out). 

Print  the  last  overshoot  stored  in  LIMITBUF  (symbol  888) 

Stop  real  time  clock. 

Print  position  in  phase  trajectory  where  it  was  stopped  by  the 
real  time  clock  (symbol  999) • 
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88.  Return  to  "beginning  to  start  new  problem  if  new  parameters 
remain  to  be  evaluated.  Unconditional  stop  if  all  parameters 
have  been  evaluated. 

89.  Has  the  required  number  of  overshoots  been  evaluated  prior  to 
starting  the  averaging  routine? 

90.  Add  the  just  evaluated  to  the  previous  sum  for  averaging. 

91.  Add  one  to  the  count  prior  to  averaging. 

92.  Clear  time  clock. 

o 

93*  Print  time,  and  ac. 

94.  Required  number  of  total  print-outs? 

95-  Clear  buffers. 

96.  Obtain  the  average  of  ©c- 
97-  Subtract  1.0. 

98.  Print  the  average  of  1.0  and  how  many  overshoots  were 

utilized  to  obtain  the  average. 

99*  Is  less  than  1.005  radians? 

100.  Clear  buffers. 

101.  Start  real  time  clock. 

102.  EXIT  (automatically  set  to  desired  jump-out). 

103.  Has  30  seconds  of  real  time  elapsed  since  the  last  print-out? 
10^4- .  Print  the  position  in  the  phase  trajectory  where  it  was  stopped 

by  the  real  time  clock  and  being  less  than  1.005  radians 

(symbol  1005 )  • 
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Jump  Switches: 


-  For  LIMITCYC ,  allows  the  maximum  overshoot  print- out  only. 
With  down,  a  phase  trajectory  is  printed  for  a  maximum 
of  12  cycles. 

-  Set  the  number  of  print-outs  depending  on  zeta.  Change 
the  number  of  overshoot  values  used  to  obtain  the  average 
value. 

J  -  Use  the  real  time  clock  for  an  automatic  recycle  to  a 

3 

new  parameter  when  have: 

a.  ll/2  minutes  maximum  after  the  last  print-out. 

b.  30  seconds  maximum  and  is  1*005  radians  or  less. 


Stop  Switches: 

S^  -  Stops  at  the  end  of  this  parameter  run  or  stops  at  the  end 
of  the  manual  print-out. 
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APPENDIX  B 


COMPUTER  PROGRAM 


00070 

ORG 

00007 

00007 

74 

0 

FAULT  EXF 

0 

00070 

75 

0 

40701 

SLJ 

0 

FAULPRIN 

ORG 

40000 

40000 

74 

0 

00100 

START  EXF 

0 

00100 

40001 

75 

0 

40004 

SLJ 

0 

MEW Z ETA- 1 

12 

1 

40623 

LA2Y2ETA  LDA 

1 

ZSTAINDX 

20 

0 

40425 

sta 

0 

ZETA 

40002 

50 

1 

00000 

ENI 

1 

0 

40003 

75 

0 

40000 

SLJ 

•  0 

START 

75 

0 

40674 

MANSTART  SLJ 

0 

MANUAL 

50' 

0 

00000 

EMI 

0 

0 

4OOO4 

12 

0 

40565 

LDA 

0 

NEWSTOP 

20 

0 

40563 

STA 

0 

STOPNOW 

40005 

12 

0 

40425 

NEWZETA  LDA 

0 

ZETA 

22 

0 

40564 

AJP 

0 

STOP  I NDX 

40006 

32 

0 

40665 

FMU 

0 

TWO 

32 

0 

40426 

FMU 

0 

OMEGAN 

40007 

20 

0 

40404 

sta 

0 

A 

12 

0 

40624 

LDA 

0 

ZETAINDX+1 

40010 

31 

0 

40425 

FSB 

0 

ZETA 

40011 

22 

0 

40433 

AJP 

0 

OUIKCOMP 

22 

2 

4Q433 

Ao? 

2 

QUIKCOMP 

12 

0 

40402 

LDA 

0 

POINT01 

40012 

20 

0 

40154 

STA 

0 

TABLES-i-l 

20 

0 

40173 

STA 

0 

TABLEL+1 

40013 

20 

0 

40207 

STA 

0 

TABLEM+1 

12 

0 

40661 

LDA 

0 

COUNT 10 

40014 

20 

0 

404  00 

STA 

0 

TENTHSEC 

50 

0 

00000 

ENI 

0 

0 

40015 

75 

2 

40437 

ZZ  SLJ 

2 

BIGZETA 

50 

0 

00000 

EMI 

0 

0 

40016 

12 

0 

40652 

40615 

LDA 

0 

STOP  20 

20 

0 

STA 

LDA 

0 

STOPZETA 

40017 

12 

0 

406o2 

0 

COUNT  12 

20 

0 

40617 

sta 

0 

COUNTZET 

40020 

12 

0 

4067? 

LDA 

0 

EIGHT 

20 

0 

40620 

STA 

0 

DIVIDE 

40021 

10 

0 

00010 

ENA 

0 

10 

61 

0 

40562 

SAL 

0 

STOPRINT 

40022 

75 

A 

71000 

GOAROUND  SLJ 

lx 

DECOF 

50 

6 

00000 

ENI 

6 

0 

40023 

02 

0 

40427 

02 

0 

FORZEROS 

04 

0 

00000 

04 

0 

0 

40024 

12 

0 

40664 

LDA 

0 

ONE 

20 

0 

40156 

STA 

0 

UDOT 

40025 

10 

0 

00000 

ENA 

0 

0 

20 

0 

00000 

STA 

0 

0 

40026 

20 

0 

40162 

STA 

0 

THETA 

20 

0 

40161 

STA 

0 

THETADOT 

40027 

20 

0 

40157 

STA 

0 

U 

20 

0 

40155 

STA 

0 

1 

40030 

20 

0 

40377 

STA 

0 

INDEX 

20 

0 

40615 

STA 

0 

' NDEXAVE 

40031 

20 

0 

40351 

STA 

0 

PRI  NTBUF-r3 

20 

0 

40352 

STA 

0 

PRINTBUF+4 

40032 

20 

0 

40353 

STA 

0 

PRI NTBUF+5 

20 

0 

40614 

STA 

0 

LIMINDEX 

40033 

20 

0 

40607 

STA 

0 

LIMITBUF+2 

20 

0 

40610 

S  l  A 

0 

LIMITBUF+3 

40034 

20 

0 

40611 

sta 

0 

LI M I TBUF+4 

20 

0 

40612 

STA 

0 

LIMITBUF+5 

40035 

20 

0 

40613- 

sta 

0 

L I M I TBUF+6 

12 

0 

40621 

LDA 

0 

PRINTSET 

40036 

20 

0 

40326 

STA 

0 

0K2PRINT+1 

20 

0 

403  35 

STA 

0 

B0UNPRIN+4 

40037 

12 

0 

40622 

LDA 

0 

NEWLIMPT 

20 

0 

40544 

STA 

0 

LIMPRINT+1 

4o040 

12 

0 

40253 

LDA 

0 

FLOTDAT  A  . 

- 

20 

0 

40071 

STA 

0 

FLOATEST 

40041 

12 

0 

40254 

LDA 

0 

FLOTDATA-i-1 

20 

0 

40072 

sta 

0 

FLOATEST+1 

40042 

75 

0 

40452 

NEWINPUT  SLJ 

0 

I NDXVALU 

50 

0 

00000 

ENT 

0 

00000 

40043 

12 

0 

404io 

LDA 

Q 

DELTA 

40044 

22 

0 

0 

40500 

AJ? 

0 

LINEAR 

12 

40513 

LDA 

0 

NOSINK 

40045 

20 

75 

0 

4 

Toll 

S  i  A 
SLJ 

0 

4 

FLOATEST 

DECOF 

50 

0 

00000 

ENI 

0 

0 

40046 

01 

0 

40404 

01 

0 

A 

06 

0 

00001 

06 

0 

1 

40047 

72 

0 

40046 

RAO 

0 

7-1 

40050 

2? 

§ 

71000 

40412 

SW 

4 

0 

DECOF 

MO INERT 

06 

0 

00000 

06 

0 

0 

40051 

72 

0 

40510 

RAO 

0 

LINEAR+10 

75 

4 

71000 

SLJ 

4 

DECOF 

40052 

01 

0 

40420 

01 

0 

KMOTCONS 

06 

0 

00000 

06 

0 

0 

40053 

10 

0 

00000 

combined 

ENA 

0 

0 

40054 

20 

0 

40160 

STA 

0 

QX 

20 

0 

40163 

STA 

0 

QY 

75 

4 

60200 

SLJ 

4 

RUNGE 

40055 

00 

0 

40153 

.  0 

0 

tables 

00 

0 

40164 

0 

0 

DERIVES 

40056 

75 

4 

60201 

SLJ 

4 

RUNGE+1 

50 

0 

00000 

ENI 

0 

0 

40057 

75 

0 

40310 

SLJ 

0 

COMB PR IN 

50 

0 

00000 

ENI 

0 

0 

40060 

12 

0 

40155 

LDA 

0 

T 

20 

0 

40174 

STA 

0 

TL 

40061 

20 

0 

40210 

STA 

0 

tm 

12 

0 

40156 

LDA 

0 

UDOT 

40062 

20 

0 

40175 

STA 

0 

VDOT 

33 

0 

40411 

FDV 

0 

RHO 

40063 

20 

0 

40211 

STA 

0 

W  DO  T 

12 

0 

40157 

LDA 

0 

U 

40064 

20 

0 

40176 

STA 

0 

V 

40065 

20 

0 

40200 

STA 

0 

THETADL 

33 

0 

40411 

FDV 

0 

RHO 

20 

0 

40212 

STA 

0 

W 

40066 

20 

0 

40214 

STA 

0 

THETADM  • 

12 

0 

40162 

LDA 

0 

THETA 

40067 

20 

0 

40201 

STA 

0 

THETAL 

33 

0 

40411 

FDV 

0 

RHO 

40070 

20 

0 

40215 

STA 

0 

THETAM 

50 
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00024 

40653 

00 

0 

OOOOO 

ST0P17 

OCT 

21 

00 

0 

00021 

40654 

00 

0 

OOOOO 

STOP25 

OCT 

31 

00 

0 

00031 

40655 

00 

0 

OOOOO 

C0UNT4 

OCT 

4 

00 

0 

00004 

40656 

00 

0 

OOOOO 

C0UNT5 

OCT 

5 

- 

00 

0 

00005 

40657 

00 

0 

OOOOO 

OCT 

5 

00 

0 

00005 

40660 

00 

0 

OOOOO 

C0UNT7 

OCT 

7 

00 

0 

00007 

40661 

00 

0 

OOOOO 

COUNUO 

OCT 

12 

00 

0 

00012 

40662 

00 

0 

OOOOO 

COUNT 12 

OCT 

14 

00 

0 

00014 

40663 

00 

0 

OOOOO 

OCT 

14 

00 

0 

00014 

40664 

20 

0 

14000 

ONE 

DEC 

loO 

00 

0 

OOOOO 

40665 

20 

0 

24000 

TWO 

DEC 

2oO 

00 

0 

OOOOO 

40666 

20 

0 

24000 

DEC 

2*0 

00 

0 

OOOOO 

40667 

20 

0 

26000 

three 

DEC 

3  o  0 

00 

0 

OOOOO 

40670 

20 

0 

26000 

DEC 

3  o  0 

00 

0 

OOOOO 

40671 

20 

0 

34000 

FOUR 

DEC 

4a  0 

00 

0 

OOOOO 

40672 

20 

0 

44000 

EIGHT 

DEC 

8  o  0 

00 

0 

OOOOO 

40673 

20 

0 

35000 

FIVE 

DEC 

5  a  0 

00 

0 

OOOOO 

40674 

50 

1 

OOOOO 

MANUAL 

ENI 

1 

0 

75 

4 

71000 

SLJ 

4 

DECOF 

40675 

01 

04 

0 

0 

40605 

01570 

oi 

04 

0 

0 

L!MJTBUF 

40676 

50 

2 

OOOOO 

ENI 

2 

0 

75 

4 

71000 

SLJ 

4 

DECOF 

40677 

01 

0 

40346 

01 

0 

PRINTBUF 

04 

0 

03651 

04 

0 

3651 

40700 

75 

0 

40563 

SLJ 

0 

STOPNOW 

50 

0 

OOOOO 

ENI 

0 

0 

40701 

50 

1 

OOOOO 

FAULPRIN 

ENI 

1 

0 

75 

4 

71000 

SLJ 

4 

DECOF 

40702 

01 

0 

40605 

01 

0 

LIMITBUF 

- 

04 

0 

01570 

04 

0 

1570 

40703 

50 

2 

OOOOO 

ENI 

2 

0 

75 

4 

71000 

SLJ 

4 

DECOF 

4  0704 

01 

0 

40346 

01 

0 

PRINTBUF 

40705 

04 

0 

01232 

04 

0 

1232 

74 

0 

42003 

EXF 

0 

42003 

76 

0 

40000 

SLS 

0 

START 

40706 

20 

0 

14000 

RESTVALU 

DEC 

IcO 

00 

0 

OOOOO 

40707 

00 

0 

OOOOO 

DEC 

0 

00 

0 

OOOOO 

40710 

20 

0 

04631 

DEC 

0  a  6 

46 

3 

14631 

40711 

20 

0 

063  14 

DEC 

0  a  8 

63 

1 

46314 

40712 

00 

0 

OOOOO 

BSS 

14 

00 

0 

OOOOO 

40726 

17 

7 

64631 

DELTAVAL 

DEC 

a  3 

46 

3 

14631 

94 


40727 

17 

7 

46314 

63 

1 

46314 

40730 

17 

7 

27534 

12 

1 

72702 

40731 

17 

7 

15075 

34 

1 

217  27 

40732 

17 

7 

05075 

34 

1 

21727 

40733 

00 

0 

00000 

00 

0 

00000 

40747 

17 

7 

46314 

63 

1 

46314 

40750 

17 

7 

56314 

40751 

63' 

1 

46314 

20 

0 

04000 

00 

0 

00000 

40752 

20 

0 

06314 

63 

1 

46314 

40753 

20 

0 

07146 

31 

4 

63146 

40754 

00 

0 

00000 

00 

0 

00000 

40774 

00 

0 

00000 

00 

0 

OQOOO 

40775 

17 

7 

5  6314 

63 

1 

46314 

40776 

17 

7 

66314 

63 

1 

46314 

40777 

20 

0 

04631 

46 

3 

l-'-63T 

41000 

20 

0 

06314 

63 

46314 

41001 

20 

6 

14000 

00 

0 

00000 

41002 

17 

7 

35075 

34 

1 

21727 

41003 

17 

7 

46314 

63 

1 

46314 

41004 

20 

5 

07146 

" 

31 

4 

63146 

41005 

00 

0 

00000 

00 

0 

00000 

DEC 

9  1 

DEC 

«  03 

DEC 

0  a  0 1 

DEC 

0*005 

BSS 

14 

JM07VALU 

DEC 

6  1 

DEC 

a  2 

DEC 

o  5 

DEC 

08 

"EC 

o9 

BSS 

20 

FLFTRATO 

DEC 

0 

DEC 

•  2 

DEC 

c4 

DEC 

0  6 

DEC 

68 

DEC 

leO 

DEC 

0o04 

DEC 

0  c  1 

DEC 

0  0  9 

BSS 

20 

END 

95 


:>olut\or\  3 

2$fc)n99  999999 

2a  99999999 

79999999 

e  99999999 

A  29999999 

(O  99999999 

vtuwi  be  r 

Tu 

Jm  49999999 

Tl  49999999 

J^Tl99999999 

(-^  20000000 

fL  79999999 

p  99999999 

K  99999999 

X  19999999 

1 

Jj»  40000000 

15999999 

1 

OJn  99999999 

5  49999999 

po'»y\-t  1 

7%  U  hkbef*  2 

10000000 

95004083- 

1 

■>  n  48334152- 

2 

T, 

1 

19999999 

*  18006400 

18669244- 

1 

■J  L 

3 

4 

5 

t  imim 
50000000 

©c  25531729 

mn% 

©t  40519239- 

\mnr 

1 

1 

6 

59999999 

42475713 

14458358 

7 

69999999 

46361850 

18907183 

8 

79999999 

49437347 

23703703 

9 

89999999 

51750106 

28769225 

10 

99999999 

53350719 

34029984 

11 

10999999 

1 

54291737 

39417381 

12 

11999999 

1 

54626982 

44868139 

13 

12999999 

1 

54410913 

50324400 

14 

14000000 

1 

53698048 

55733762 

15 

15000000 

1 

52542443 

61049253 

16 

16000000 

1 

50997217 

66229265 

17 

17000000 

1 

49114140 

71237437 

18 

18000000 

1 

46943268 

76042502 

19 

19000000 

1 

44532629 

80618101 

20 

20000000 

1 

41927962 

84942563 

21 

21000000 

1 

39172497 

88998670 

22 

22000000 

1 

36306785 

23 

23000000 

1 

33368563 

24 

24000000 

1 

30392671 

99445851 

25 

25000000 

1 

27410989 

10233595 

1 

26 

26000000 

1 

24452419 

10492882 

1 

27 

27000000 

1 

21542892 

10722807 

1 

28 

28000000 

1 

18705407 

10923980 

1 

• 

29 

29000000 

1 

15989086 

11097264 

1 

0™  15926240 

0™  11097173 

1 

30 

30000000 

1 

13625000 

11245019 

1 

12991268 

11241908 

1 

31 

31000000 

1 

11610459 

98937811- 

11370928 

1 

99052498- 

1 

11356496 

1 

32 

32000000 

1 

1 

11478221 

1 

67134562- 

1 

11439661 

1 

33 

33000000 

1 

84309241- 

1 

11569649 

1 

34535657- 

1 

11490539 

1 

34 

34000000 

1 

71843596- 

1 

11647559 

1 

15682165- 

2 

11508610 

1 

35 

35000000 

1 

61221074- 

1 

11713950 

1 

-31509827- 

1 

11493639 

1 

36 

36000000 

1 

52169158- 

1 

11770525 

1 

-64486655- 

1 

11445624 

1 

37 

37000000 

1 

44455624- 

1 

11818734 

1 

-97189429- 

1 

11364757 

1 

38 

38000000 

1 

37882584- 

1 

11859816 

1 

-12947823 

11251384 

1 

39 

39000000 

1 

32281408- 

1 

11894823 

1 

-16124091 

11105976 

1 

40 

40000000 

1 

27508401- 

1 

11924654 

1 

-19238868 

10929107 

10721427 

1 

41 

41000000 

1 

23441113- 

1 

11950075 

1 

-22285244 

1 

42 

42000000 

1 

19975199- 

1 

11971737 

1 

-25257956 

10483648 

1 

43 

43000000 

44000000 

1 

17021742- 

1 

11990196 

1 

-28153122 

10216527 

1 

44 

1 

14504971- 

1 

12005926 

1 

-30968015 

-33700866 

99208537 

45 

45000000 

1 

12360321- 

1 

12019330 

1 

95974404 

46 

46000000 

1 

10532771- 

1 

12030752 

1 

-36350703 

92471132 

C  or,  tact  1 

46600000 

1 

95686438- 

2 

12036778 

1 

-37900604 

90243442 

yvvi  m  otr  1 

46600000 

1 

-37900604 

12036778 

1 

95686438- 

2 

90243442 

47 

47000000 

1 

-35550956 

11889925 

1 

-60605683- 

2 

90374327 

48 

48000000 

1 

-30294526 

11561398 

1 

-39333994- 

1 

90140801 

89611697 

49 

49000000 

1 

-25815292 

11281446 

1 

-65385515- 

1 

50 

50000000 

1 

-21998341 

11042886 

1 

-85373205- 

1 

88853284 

51 

51000000 

1 

-18745749 

10839599 

1 

-10028037 

87921151 

52 

52000000 

1 

-15974074 

10666370 

10518753 

1 

-11094153 

86861815 

53 

53000000 

1 

-13612208 

1 

-11806455 

85714102 

APPENDIX  C 

PHASE  TRAJECTORY  PRINT  OUT 

(Numbers  indicated  in  tenths  followed  by  power  of  10 ) 


Salutio^  1 

ZSlO,J 9999999 

3*^11111111 

K/fT  39999999 

ftu  w\  ler 

Jm  99999999- 

1 

Tu  89999999 

Twit  pi1 99999999 

K  99999999 

)C /Tm99999999 

1 

Pm/  Tm47999999 

1 

humber1 

IT  33800000 

1 

ec  52456182- 

1 

©c  13396964 

1 

of  2 

11230000 

2 

38509278- 

1 

11053498 

1 

3 

21289999 

2 

12863485- 

1 

10697107 

1 

4 

32499999 

2 

79971808- 

2 

10583527 

1 

5 

44279999 

2 

65238396- 

2 

10542905 

1 

6 

56369999 

2 

59130680- 

2 

10524201 

1 

7 

68649998 

2 

55271237- 

2 

10511522 

1 

8 

81049998 

2 

53208045- 

2 

10504778 

1 

9 

93519998 

2 

51991503- 

2 

10500457 

1 

10 

10602999 

3 

51279129- 

2 

10498432 

1 

3  numlotro4mOTt>v7\uw>$ 

a  Vera  o,  e  c\ 


e  99999999 

A  29999999  p  99999999 

fw  47999999 

31999999  p"1  99999999 

fu/  JL  35555555 

99999999  ^  39999999 

50122288- 

1  flytras.c  Vosit'ivt 
^  rrar  ^  n  L  i  yyv  t  t  C.  ^  c  1 

APPENDIX  D  1 

TYPICAL  PRINT  OUT  ILLUSTRATING  LIMIT  CYCLE 
(Numbers  indicated  in  tenths  followed  "by  power  of  10) 


2 

2$Wn 79999999 

TU/T,.  11111111 

fL//=T  59999999 

C  99999999 

A  29999999 

P  99999999 

3»w  99999999- 

1 

Tl  89999999 

99999999 

32000000 

■Fu  47999999 

p1 99999999 

K  99999999 

fc/J*,99999999 

1 

{m/Tr*  32000000 

1 

fi_ J  Ji  53333333 

U)^  99999999 

39999999 

1 

-t  35400000 

1 

et  17098208- 

1 

©c.  13104190 

1 

0<*  7 

11580000 

2 

17034145- 

1 

10721242 

1 

3 

23279999 

2 

18292628- 

2 

10323347 

1 

4 

38859999 

2 

35712402- 

3 

10200613 

1 

5 

58729999 

2 

72447348“ 

4 

10138780 

1 

6 

83939998 

2 

11243587- 

4 

10098946 

1 

7 

11619999 

3 

10685365- 

5 

10071240 

1 

8 

15810000 

3 

47854719- 

7 

10051391 

1 

9 

21077001 

3 

29068282- 

8 

10037079 

1 

1005 

215160011794  3 

-167815371770  - 

2 

r 

999996862854 

100  5*  \r\c\\C(\it-$  6-cmax  less  tUjiA  \  QOS'  0 >r\i  zo  sec.  ot 


realty  we.  W  e|Q(>6<4  siv\ec  tv.e  'cict  A^x’.muvn 
UJas  Hau)  s e\ot'~y } • 

TKe-  t;^vw?  Velocity,  DAcl  ^  tVie  i  Cl  it  Cow^t^d 

*f>o\vVt  i*  4. 


APPENDIX  D  2 

TYPICAL  PRINT  OUT  ILLUSTRATING  NO  LIMIT  CYCLE 
(Numbers  Indicated  in  tenths  followed  hy  power  of  10 ) 


